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*n , Abstract. We present a detailed study of chemical freeze-out in nucleus- 

nucleus collisions at beam energies of 11.6, 30, 40, 80 and 158^4 GeV. By 
, analyzing hadronic multiplicities within the statistical hadronization approach, 

we have studied the chemical equilibration of the system as a function of cen- 
ter of mass energy and of the parameters of the source. Additionally, we have 
tested and compared different versions of the statistical model, with special 
emphasis on possible explanations of the observed strangeness hadronic phase 

2 ■ space under-saturation. 

Ch ■ 
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1. The Statistical Hadronization Model 

One of the main results of the study of high energy A-A collisions is a surprising 
success of the statistical-thermal models in reproducing essential features of particle 
production (see for example [HE]). In this paper we study nucleus- nucleus collision 
within the statistical model in the energy range that is believed [ |3] to cover the 
threshold for creation of Quark-Gluon Plasma (QGP) in the early stage of Pb-Pb 
collisions. 

The main idea of the SHM is that hadrons are emitted from regions at statistical 
equilibrium. No hypothesis is made about how statistical equilibrium is achieved; 
this can be a direct consequence of the hadronization process. In a single collision 
event, there might be several clusters with different collective momenta, different 
overall charges and volumes. However, Lorentz-invariant quantities like particle 
multiplicities are independent of clusters momenta. 
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If the system size is sufficiently large [U], the analysis can be done in grand- 
canonical ensemble in which the mean primary multiplicity of the j th hadron with 
mass rrij and spin Jj reads: 



where T is the temperature, V the scaling volume, = (Qj, Bj, Sj) is a vector 
having as components the electric charge, baryon number and strangeness of the 
hadron and fi = ()J.q, hb, Ms) is a vector of the corresponding chemical potentials; 
the upper sign applies to fermions, the lower to bosons. 

In order to correctly reproduce the data, it is also necessary to introduce at 
least one non-equilibrium parameter suppressing hadrons containing valence strange 
quarks, 75 7^ 1 [E]- With this supplementary parameter, hadron multiplicity is as 
in Eq. (JJJ with the replacement: exp[/i • q_j/T] — * exp[/i, • q.,-/T]7g s where n s stands 
for the number of valence strange quarks and anti-quarks in the hadron j. 

Finally, the overall multiplicity to be compared with the data, is calculated as 
the sum of primary multiplicity JQ) and the contribution from the decay of heavier 
hadrons: (rij) — (ni) primary + J\ Br(fc — * j)(nk), where the branching ratios are 
taken from the latest issue of the Review of Particle Physics [ Ej • 

2. Experimental Data Set And Analysis Results 

The bulk of the experimental data consists of measurements made by NA49 collab- 
oration in central Pb-Pb collisions at beam momenta of 30, 40, 80 and 158A GeV [ 
□ El. As far as AGS data at U.6A GeV [03 QUI EH HI is concerned, we have used 
both multiplicities measured by the experiments and extrapolations of measured 
rapidity distributions (for details, see [I15|). 

The analysis has been carried out by looking for the minima of the 

/ exp ^thco%2 

X 2 = Yli — ~^ — • The fitted parameters within the main scheme SHM(7s) 
are shown in tabled The observed differences in the fit parameters between two 
independent analyses A and B are of the order of the fit errors. 

A major result of these fits is that 7s is significantly smaller than 1 in almost all 
cases, that is strangeness seems to be under-saturated with respect to a completely 
chemically equilibrated hadron gas. This confirms previous findings [ 1121 1131 114| . 

2.1. Strangeness correlation volume 

To account for the observed under-saturation of strangeness, a picture has been put 
forward in which strangeness is supposed to be exactly vanishing over distances less 
than those implied by the overall volume V [EL i- e - all clusters or fireballs emerge 
with S = and they are not allowed to share non-vanishing net strangeness. 

Assuming, for sake of simplicity, that all clusters have the same typical volume 
V c and that treating baryon number and electric charge (but not strangeness) grand- 
canonically, gives a good description of the system, one can perform the analysis in 
strangeness canonical ensemble (for details, see [I15|). 

If V c is sufficiently small, the multiplicities of strange hadrons turn out to be 
significantly suppressed with respect to the corresponding grand-canonical ones due 
to an effect called canonical suppression. 
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We have fitted the data sample of full phase space multiplicities in Pb-Pb 
collisions at 158^4 GeV fixing 75 = 1. The quality of the fit is worse (x 2 /dof = 
37.2/9) with respect to the SHM(7s) model, whilst thermal parameters T=157.9 
MeV and /^s=261.5 MeV are compatible with the main version of the statistical 
model. Our result suggests that, for the local strangeness correlation to be an 
effective mechanism, the cluster volume should be of the order of 2.5% of the overall 
volume. 

2.2. Superposition of NN collisions with a equilibrated fireball 

In this picture, henceforth referred to as SHM(TC), the observed hadron production 
is approximately the superposition of two components: one originated from a large 
fireball at complete chemical equilibrium at freeze-out, with 75 = 1, and another 
component from single nucleon-nucleon collisions. Since it is known that in NN 
collisions strangeness is strongly suppressed [Q], the idea is to ascribe the observed 
under-saturation of strangeness in heavy ion collisions to the NN component. 

With the simplifying assumption of disregarding subsequent inelastic collisions 
of particles produced in those primary NN collisions, the overall hadron multiplicity 
can be written then as (nj) = (N c )(rij)]\rN + ( n j)v, where (iij)jviv is t ne average 
multiplicity of the j th hadron in a single NN collision, (N c ) is the mean number of 
single NN collisions and (rij)v is the average multiplicity of hadrons emitted from 
the equilibrated fireball. 

To calculate (rij)j\rjv we have used the statistical model and fitted pp full phase 
space multiplicities measured at the same beam energy by the same NA49 experi- 
ment. For np and nn collisions, the parameters of the statistical model determined in 
pp are retained and the initial quantum numbers are changed accordingly. Theoret- 
ical multiplicities have been calculated in the canonical ensemble, which is described 
in detail in ref. 

We have fitted T, V, /j,b of the central fireball and (N c ) by using NA49 data in 
Pb-Pb collisions at 158^4 GeV. The fit quality, as well as the obtained values of T 
and /ib, are comparable to the main fit within the SHM(7s) model. The predicted 
number of "single" NN collisions is about 50 with a 16% uncertainty. 

2.3. Non-equilibrium of hadrons with light quarks 

In this model [ 117) two non-equilibrium parameters are introduced for the different 
types of quarks, j„ for u, d quarks and j s for strange quarks. By defining: 75 = ^ 

)q 

and V = V~fq, the Boltzmann limit of average multiplicity reads: 



where Bj is the baryon number, as long as mesons have two and baryons have three 
valence quarks. The parameter 7 ? has a definite physical bound for bosons which can 

be obtained by requiring the convergence of the series 5Ziv =0 (7q * ) ex P(~ Ne/T + 
N/j, ■ cjj/T) for any value of the energy. If 7 9 reaches its bounding value 
j q = exp(m T o /2T) ~ 1.5 for T ~ 160 MeV, a Bose condensation of particles in the 
lowest momentum state sets in. 

It is seen that the absolute \ 2 minimum falls in the region of pion condensation, 
at j q ~ 1.62, with x 2 — 13 and T ~ 140 MeV. This finding is in agreement with 
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what is found in ref. [ I17j. However, there is also a local minimum at the lower 
edge 7 9 = 0.6, with a temperature of 187 MeV, which is only one unit of % 2 higher 
than the absolute minimum. This indicates that the absolute minimum could be 
rather unstable against variations of the input data and this is in fact what we find 
by varying down the pion multiplicities by only 1 a. For this case, the absolute 
minimum of \ 2 now lies at j q = 0.6 instead of at the pion condensation point. 

In view of the instability of the fit, and of the small relative x 2 improvement in 
comparison with the main fit, we conclude that there is so far no evidence for the 
need of this further non-equilibrium parameter. 

3. Energy Dependence 

The chemical freeze-out points in the /ib — T plane are shown in fig. ^ The four 
points at beam energies of 11.6, 40, 80 and 158A GeV have been fitted with a 
parabola: T — 0.167 — O.ISS/j,^, where T and [is are in GeV. 

A possible indication of deconfinement phase transition in Pb-Pb collisions at 
the low SPS energies was reported on the basis of the observed energy dependence 
of several observables [ 118) . Particularly, the (K + )/(ir + ) ratio shows a peaked 
maximum at about 30yl GeV. One may expect that this anomaly should be reflected 
in the energy dependence of 7s parameter fitted within SHM(7s) scheme. This 
dependence is plotted in fig. [21and in fact a maximum shows up at 30^4 GeV. 

The anomalous increase of relative strangeness production at 30A GeV can be 
seen also in the Wroblewski variable As = 2(ss)/((uu) + (dd)), the estimated ratio 
of newly produced strange quarks to u, d quarks at primary hadron level, shown in 

fig-El * 

In order to further study strangeness production features, we have also com- 
pared the the measured (K + )/(ir + ) ratio with the theoretical values in a hadron 
gas along the interpolated freeze-out curve for different values of 7s (see fig. 0. 
The calculated dependence of (K + ) / (tt + ) on [Ib is non-monotonic with a broad 
maximum at /is ~ 400 MeV (i.e. Ef, eam ~ 30 A GeV). Taking into account that 
systematic errors at different energies in Pb-Pb collisions are fully correlated, we 
can conclude that the data points seem not to follow the constant 7s lines. 

4. Summary and Conclusions 

It is found that the main version of the statistical model (with 7s), fits all the data 
analyzed in this paper. We have tested a model (SHM(TC)) which can fit the data 
at 158A GeV very well if the number of independent NN collisions is around 50 
with a sizeable uncertainty. A model in which strangeness is assumed to vanish 
locally [E] yields a worse fit to the data with respect to SHM(7s) and SHM(TC). 
Moreover, we have found that the present set of available data does not allow to 
establish whether a further non-equilibrium parameter (7^) is indeed needed to 
account for the observed hadron production pattern. 

The evolution of the freeze-out temperature and baryon-chemical potential is 
found to be smooth in the AGS-SPS-RHIC energy range. (K+) / (ir + ) ratio, Wrob- 
lewski factor As as well as 7s parameter calculated within the statistical model 
suggests that there might be a peak in relative strangeness production at about 
30A GeV of beam momenta. 
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Parameters 


Main analysis A 


Main analysis B 




Main analysis A 


Main analysis B 




Au-Au 11. 6 A Gev 




Pb-Pb 30A GeV 


T (MeV) 
Ms (MeV) 

7S 
V 


118.1±3.5 (4.1) 
555±12 (13) 
0.652±0.069 (0.079) 
1.94±0.21 (0.24) 


119.1±4.0 (5.4) 
578±15 (21) 
0.763±0.086 (0.12) 
1.487±0.18 (0.25) 




139.5 
428.6 
0.938±0.078 (0.13) 
6.03±0.50 (0.85) 


140.3 

428.7 
1.051±0.103 (0.16) 
5.273±0.526 (0.80) 


X^/dof 


4.0/3 


5.5/3 




5.75/2 


4.6/2 




Pb-Pb 40A GeV 




Pb-Pb 80A GeV 


T (MeV) 
Ms (MeV) 

7S 

V 


147.6±2.1 (4.0) 
380.3±6.5 (13) 
0.757±0.024 (0.046) 
8.99±0.37 (0.71) 


145.5±1.9 (3.5) 
375.4±6.4 (12) 
0.807±0.025 (0.047) 
8.02±0.34 (0.63) 




153.7±2.8 (4.7) 
297.7±5.9 (9.8) 
0.730±0.021 (0.035) 
15.38±0.61 (1.0) 


151.9±3.4 (5.4) 
288.9±6.8 (11) 
0.766±0.026 (0.042) 
14.12±0.65 (1.1) 


X^/dof 


14.7/4 


13.6/4 




11.0/4 


10.4/4 




Pb-Pb 158A GeV 








T (MeV) 
Ms (MeV) 

7S 

V 


157.8±1.4 (1.9) 
247.3±5.2 (7.2) 
0.843±0.024 (0.033) 
21.13±0.80 (1.1) 


154.8±1.4 (2.1) 
244.5±5.0 (7.8) 
0.938±0.027 (0.042) 
18.46±0.69 (1.1) 








X^/dof 


16.9/9 


21.6/9 









Table 1. Summary of fitted parameters (V'=VT 3 exp[— 0.7 GeV/T]) at AGS and 
SPS energies in the framework of the SHM(7s) model. The re-scaled errors (see [E0 
I15P are quoted within brackets. For Pb-Pb at 30 A GeV data, we have constrained 
T and ijlb to lie on the fitted chemical freeze-out curve (see fig. 0) 
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Pb-Pb SPS Au-Au AGS 




Chemical freeze-out 
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Fig. 1. LEFT: Chemical freeze-out points in the [Ib — T plane in various heavy 
ion collisions. The full round dots refer to Au-Au at 11.6 and Pb-Pb collisions at 
40, 80, 158A GeV obtained in the analysis A, whilst the hollow square dot has been 
obtained in ref. [ I19| by using particle ratios measured at midrapidity in Au-Au 
collisions at \fs NN = 130 GeV. The hollow round dot without error bars refers to 
Pb-Pb collisions at 30A GeV and has been obtained by forcing T and /is to lie on 
the parabola fitted to the full round dots. 

RIGHT: Measured (K + ) / (tt + ) ratio as a function of the fitted baryon-chemical po- 
tential. The full square dot is a preliminary full phase space measurement in Au-Au 
collisions at \/~s NN — 200 GeV [GO]- For the SPS energy points the statistical er- 
rors arc indicated with solid lines, while the contribution of the common systematic 
error is shown as a dotted line. Also shown the theoretical values for a hadron gas 
along the fitted chemical freeze-out curve (left), for different values of 75. 
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Fig. 2. LEFT: Strangeness non-equilibrium parameter 75 as a function of the 
nucleon-nucleon centre-of-mass energy. Full dots refer to fit A, hollow dots fit B. 
RIGHT: As estimated as a function of the fitted baryon-chemical potential. Also 
shown the theoretical values for a hadron gas along the fitted chemical freeze-out 
curve shown in fig-Q for different values of 75. 



